Helical instability of charged vortices in layered superconductors 
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It is shown that the electric charge of vortices can result in a helical instability of straight vortex 
lines in layered superconductors, particularly Bi-based cuprates or organic superconductors. This 
instability may result in a phase transition to a uniformly twisted vortex state, which could be de- 
tected by torque magnetometry, neutron diffraction, electromagnetic or calorimetric measurements. 
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Vortices in superconductors carry the quantized mag- 
netic flux <fio = 2 x 10 -7 Oe-cm 2 resulting from the macro- 
scopic phase coherence of superconducting state. Vor- 
tices also carry a non-quantized electric charge q caused 
by the suppression of superconductivity in the vortex core 
[lH6|. In low-T c s-wave superconductors this charge is 
usually negligible and does not manifest itself in the elec- 
tromagnetic response of vortices driven by the Lorentz 
force of superconducting currents. However, the situ- 
ation changes in superconductors with short coherence 
length £, low superfluid density and unconventional pair- 
ing symmetry combined with the competition of super- 
conductivity with non-superconducting spin or charged 
ordered states, as characteristic of high-T c cuprates, re- 
cently discovered oxypnictides or organic superconduc- 
tors pj. For cuprates, theoretical estimates [l|-|4( predict 
a relatively large fraction ~ 1CP 3 of the electron charge 
e per each pancake vortex residing on the ab planes, yet 
even larger charge of different sign was observed by nu- 
clear quadrupole resonance [|| . It has been suggested [l| 
that the vortex charge could change the sign of the Hall 
coefficient observed in cuprates Q or result in structural 
transformations of the vortex lattice [l(| ■ 

In this paper we show that vortex charge can cause 
an intrinsic helical instability of a rectilinear vortex and 
a phase transition to a twisted vortex state. This in- 
stability is different from the helical instability of vor- 
tices driven by either currents flowing along the vortex 
line or by screw dislocations [l3| or twisted vortex 
states in rotating liquid He 14|. The buckling insta- 
bility of vortices results from the Coulomb repulsion of 
charged pancake vortices which tend to shift away from 
the straight line along the c-axis as illustrated by Fig. 1. 
Such charge fragmentation is inhibited by the vortex line 
tension caused by weak magnetic and Josephson coupling 
of vortex pancakes 17, 3], and also by charge screening, 
which confines the relative displacements of pancakes on 
neighboring ab planes within the Thomas-Fermi screen- 
ing length Xtf- Thus, the helical instability would be 
most pronounced in layered materials with low vortex 
line tension and Xtf ~ as characteristic of high-T c 
cuprates, ferropnictides or organic superconductors. 

To calculate properties of spiral vortices we write the 





FIG. 1: Spiral instability of a straight chain of charged pan- 
cake vortices in a layered superconductor. 



excess linear charge p(r) in a vortex as follows 



p(r) 



+ Pa exp(-r 2 /2£ 2 ) 



(1) 



Here the first term is the BCS contribution resulting from 
the change in the chemical potential p around the core, 
p BCS {r) oc [A 2 (r) - A 2 ], A(r) s A r/(r 2 + £ 2 ) 1/2 is the 
modulus of the order parameter, po = eNAqdhiT c /dp, 
N is the density of states at the Fermi surface in 
the normal state. The BCS vortex charge qbcs — 
27r£ 2 po hi(A/£) is spread over the London penetration 
depth A [2j|. Strong dependence of the critical temper- 
ature T c on doping enhances po in cuprates. The term 
oc p a in Eq. (pj is added phenomenologically to take 
into account the localized core charge due to competing 
superconducting and antiferromagnetic orders in uncon- 
ventional superconductors [!, Q . NMR experiments indi- 
cate [|| that the local core charge in cuprates can greatly 
exceed the BCS contribution. Eq. ([T]) corresponds to 
the following Fourier transform p(k) = 2ir£ i 2 [poKo(k£ i ) + 
p a exp(— fc 2 £ 2 /2)] and the total excess charge per unit 
length q ~ 27r£ 2 [p ln(A/£) + p a ] where K (x) is the mod- 
ified Bcssel function. 

The excess charge density p{r) in a superconductor is 
screened in the same way as in a normal metal [3-3]. 
Screening is determined by the Fourier transform of the 
static dielectric function e(k), which, for the isotropic 
Thomas- Fermi model equals e(k) = 1 + K 2 /k 2 where 
k = 1 /Xtf- The Fourier transform of the electric poten- 
tial (p(k,p) produced by a charged curved vortex param- 
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eterized by its displacement u(z) relative to the z— axis 
is given by the Poisson equation: 



(k 2 +p 2 )e(k,p)ip(k,p) =4 7 r ( o(k) 



where e(k,p) in a uniaxial material depends on both the 
in-plane wave vector k and the z-component p perpendic- 
ular to the layers. From Eq. we obtain the functional 
of electrostatic energy VF{u(z)} = (1/2) J pipd 3 r: 



W 



\p{k)\ d \s.dpdz\dz 2 
47r 2 (fc 2 + p 2 )e(k,p) ' 



,ip(z 2 -zi)+'ik[u(zi)-u(z 2 )] 



(3) 

Here two periodic structures u(z) are considered: heli- 



cal distortions, 



u cos Qz and u y = u sin Qz, and 



planar zig-zag distortions, u x — u cos Qz and u y = 
where u and Q quantify the amplitude and the period of 
the structures. For the spiral vortex, we have ku(zi) — 
ku(z2) = 2usm{Qz-)[k x sin(Qz + ) + k y cos(Qz + )] where 
z± = (zi ± z-2)/2. Neglecting a possible dependence of 
p(k,p) on p due to charge modulation along the z-axis, 
integrating Eq. ([3]) over z\ + Z2, and the polar angle in 
the k plane, and adding the elastic energy F e , gives the 
total line energy of a vortex helix F s = F e + W s where: 

w f°°u\ (VU^r* J e- 2lpz Jo[2uk\sinQz\} 
W s = J o k\ P (k)\ dkjjpdz 77{k2+p2)e{kjp) , 



„ u 2 e Q 2 



2„,2 



A 2 7 



£ 2 (1 + A 2 Q 2 ) 4A 



(4) 

+ ^ln(l + A 2 Q 2 ) (5) 



Here F e describes the dispersive tilt energy of a vortex in 

Jq(x) is the Bessel 



17, m 



a uniaxial superconductor 
function, 7 = A c /A is the anisotropy parameter, and 
£0 = ( ( / ) o/47rA) 2 is the vortex energy scale. For a zig- 
zag vortex, we obtain F z = F e /2 + W z , where W z is 
given by Eq. (|4]) in which Jo[2uk\ smQz\] is replaced by 
Jq [uk\ sinQz|]. To determine which of the two structures 
has lower energy, we minimize F s and F z with respect to 
u and Q using e(k,p) for a layered metal [lij : 



e(k,p) = e + 



60S k sinh(fcs) 
2[cosh(fcs) — cos(ps)]sk ' 



k < 2k F (6) 



where k 2 = 4n 2 e 2 h 2 / m* e^s 2 , s is the interlayer spacing, 
6q is the background dielectric constant, m* is the elec- 
tron effective mass, and kp is the Fermi momentum. For 
k > 2k p, the last factor k in the denominator should be 
replaced by k — ^Jk 2 — Akp. Eq. ([6]) takes into account 
the anisotropy of screening at large fc, and the Friedel 
oscillations due to singularity in de/dk at k = 2kp. For 
(ks,ps) <C 1, Eq. ([6]) reduces to the Thomas-Fermi di- 
electric function e(k,p) = [1 + n 2 /(k 2 + p 2 )]eo with the 
screening length Xtf = k 1 - 

Now we show that a rectilinear vortex along the c- 
axis becomes unstable with respect to bending distortions 
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FIG. 2: Energy of the helical vortex line as a function of the 
wave vector Q and different ratios of q/q c described by eq. 
© for A/s = 10 3 and (2A T f/s) 2 = 10. 



if q exceeds a critical line charge q c . At the instability 
threshold q sa q c , Eq. (|H) can be expanded in small u, 
and the z— integration produces the S functions at k = 
and k = ±Q yielding the following change in W e 



5W S 



u 

T 



\ P (k)\ 2 k 3 



\p(k)\ 2 k 

e(fc,0) {k 2 + Q 2 )e(k,Q) 



dk (7) 



Hence helical distortions do reduce W e , the electrostatic 
energy gain increasing as Q increases. The quadratic ex- 
pansion of Eq. (0J for zig-zag distortions yields SW Z = 
SW s /2. Given that the charged vortex core is typi- 
cally larger than either Xtf and s, we expand e(k,p) 
in (ks) 2 <C 1 since the integral in Eq. (J7J is mostly 
determined by the region k 2 <C k 2 , and |p(fc)| 2 rapidly 
decreases for k > As the result, the energy change 
for small u takes the form 



SF„ = F e - 



q 2 u 2 



4r4 



1 - 



4sin 2 (Qs/2) 



(2Q/K) 2 sin 2 (Qs/2) 



(8) 

Here the effective core radius r§ is defined by q 2 j2r\ = 
r o °° k 3 \p(k)\ 2 dk = 4vr 2 f™(Vp) 2 rdr and Eq. Q, giving 
r o = £ a t Pa 3> Po- As q exceeds q c , the function 6F S (Q) 
shown in Fig. 2 first becomes negative at a finite Q. 
Such behavior reflects the effect of crystalline anisotropy, 
which strongly reduces the vortex line tension at QX 3> 1 
thus facilitating the short wavelength instability. The 
equation 8qSF s — at QA 3> 1 yields: 



2 = K 2 /21n( 7 /CQ), 



(9) 



so that the twist pitch i s ~ 2 3 / 2 tt \u 1/2 (j/S,k)X tf ~ 
IOAtf- For Xtf = 0.5 — 1 nm in cuprates [llj, t s ~ 5 — 10 
nm turns out to be larger than £. From the equation 
5F(q c ,Q) = and Eqs. (JH])- (EJ) we obtain the critical 
charge q c strongly reduced by crystalline anisotropy: 



2(roK) 4 eo eo[ln(7/£Q) + l/2]/7 2 



(10) 
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Given the relation 8F S = 25 F z , both helical and zig-zag 
instabilities occur at the same q c and Q, so to see which 
of these structures has lower energy, the amplitude of 
spontaneous distortions u at q > q c is to be calculated. 
Near the instability threshold q w q c , the general Eq. (UJ) 
can be expanded in powers of small u up to terms ~ u A 
and integrated at (tok) 2 ^ 1 as before. This gives the 
energy change for the spiral vortex: SF s /Fq = —a s u 2 + 
/3 s u 4 /4, where a. = (1 - q 2 Jq 2 ), & - 6Q 2 /r 2 ( K 2 + 4Q 2 ), 
and F = q 2 Q 2 /4r 2 K 2 (n 2 + Q 2 )e . Minimization of SF 
yields the dependence u(q) characteristic of the second 
order phase transition: 



2 / , K " 
3 V 4+ ^ 



(r 2 



Kr 2 



hi — 



(11) 



where C = r 2 J °° \p{k)\ 2 k 5 dk/ \p(k)\ 2 k 3 dk -> 1 if 
Pa 3> Po- For <7 ~ 9c the amplitude of the vortex 
helix is of the order of £, and the total energy gain 
equals 6F S — —a 2 F /f3 s . For a zig-zag vortex, we ob- 
tain a 2 = q s /2 and fl, = 3/3 s /8. Thus, 5F Z = 2£F s /3, 
so a helical vortex, which provides the maximum spacing 
between charged vortex pancakes at a given u, is more 
energetically favorable than a zig-zag vortex, which can 
lower its energy by transverse buckling distortions. 

The instability criterion q > q c depends on T. For ex- 
ample, the BCS vortex charge q oc Aq£ 2 in Eq. (fTJ is 
independent of T at T c — T <C T c , while q c oc £>/ 2 £ 2 c>c 
(1 -T/T c )- 1 / 2 in Eq. {TDJ diverges at T c , suggesting that 
the helical instability occurs below a certain temperature 
Th < T c . However, the NMR experiments show that 
the observed q is mostly determined by the non-BCS core 
contribution modeled by the term cx p a in Eq. (TT]). Cur- 
rently little is known about p a (T), so we analyze the 
criterion q > q c at low T where it can be expressed in 
terms of observable parameters. It is convenient to re- 
write q > q c in the form rj > rj c where rj c = qs/e is the 
fraction of the electron charge e per pancake vortex, and 



s(r K 



2 3 / 2 A7 



In 



1/2 



(12) 



For YBCO with e = 25, A TF = 0.5 nm 11], A = 200 
nm, ro = 1.5 nm, s = 0.85 nm, and 7 = 5, Eq. (|12[) gives 
rjc ~ 1.6, much larger than r\ ~ (0.2 — 2) x 10~ 2 observed 
for the optimally doped YBa2Cu3C>7 [§]. Larger values of 
7] ~ (1-5) x 10" 2 were observed for YBa 2 Cu 3 8 0. The 
situation becomes more interesting for layered cuprates 
and organic superconductors, for which 7 ~ 100 — 600 
0. For Bi-2212 with s = 1.5nm, A = 200nm, 7 = 500, 
eo = 10, and nro = 3, we obtain r\ c ~ 4 x 10~ 2 per double 
CuO planes. Therefore, layered cuprates (particularly 
underdopcd ones) and organic superconductors would be 
promising candidates for the experimental search for he- 
lical vortices, particularly at low T where the vortex core 



size in the clean limit tq(T) ~ £(0)T/T c may decrease 
due to the Kramer-Pesch effect [l6( . Such core shrinkage 
strongly reduces q c in Eq. (fT0|) and could result in an 
unusual case of ro(T) < Xtf for which the instability 
is further enhanced by stronger Coulomb interaction of 
pancake vortices. 

The single vortex helical instability may result in a 
long-range twist of the interacting vortex lattice. Indeed, 
if helical displacements u(z) of all vortices are phase 
locked, they do not change the flux density, V • u = 
and thus contribute to neither the shear nor the com- 
pression energy of the twisted vortex lattice. Thus, 
as far as the elastic and electrostatic energies are con- 
cerned, vortex structures with a long-range chiral or- 
der would be more energetically favorable than struc- 
tures with different signs of Q or phases of helical dis- 
tortions on neighboring vortices. In this case the vortex 
lattice would undergo a phase transition at q > q(T, B) 
to a uniformly twisted state. Fluctuations and pinning 
of vortices and proliferation of topological defects may 
destroy the long range chiral order at higher T and 
B, however if the spacing between pinning centers is 
much greater than the twist pitch £, pinning does not 
affect the single-vortex helical instability. The mean- 
field phase transition at q — q c results in the specific 
heat jump AC* = 2F TB(d T <x j ) 2 F / ( /> /3 s = TB{k 2 + 
4Q 2 )(d T g-d T g c )730o>« 2 (K 2 + Q 2 )e o . If d T q c » cVz, we 
have AC/ACo ~ TB^ 2 \n(-f/^K)/T c B c2 X^ F j 2 e where 
ACo = T c {8tH c ) 2 /8tt is the specific heat jump at T c . 

Interaction of vortices can be taken into account 
by adding the elastic twist energy £?(/> u 2 Q 2 /167r(l + 
Q 2 A 2 ) ~ o Bu 2 /167rA 2 [H E3 in Eq. ©. Then 
the problem reduces to the helical instability of a sin- 
gle vortex with a field-dependent line tension £; = 
(eo/7 2 ) ln(7/£Q) + (f) B/8ir\ 2 Q 2 where the last term re- 
sults from the magnetic cage potential [13, 18 1. Mini- 
mization of F(Q) at q = q c and Qs < 1 yields 

^ta^-lWtf + ^V, (13) 



9 = = 2£oe #(^ + 0=)(^lni + |^). (14) 

For B < o /87tA| f 7 2 , Eqs. {T5J} and (0 reduce 
to Eqs. (J9]) and flTDJ}. For B > (/) /87rA| F 7 2 , we 
have Q - (B/fo) 1 / 4 ^^ 1 / 2 , which gives the critical 
charge q c ~ (tq/ XXtf)(^o4'oB/4:'k) 1 ^ 2 independent of 
anisotropy. The instability region T < Th(B) defined 
by q{Th) > q c {Th,B) thus widens as B decreases. 

Helical distortions with Q ~ k can produce minibands 
in the spectrum of core quasiparticles moving along the 
vortex. This may affect the vortex viscosity, vortex mass, 
and pinning, and also smear the discrete core levels in 
the STM images of a helical vortex. Vortex chirality 
also manifests itself in a "fountain" -like currents along 
the z-axis 13J, |14|, and features of flux dynamics con- 
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FIG. 3: Mechanism of the torque exerted by a transverse 
current on the vortex helix where closed lines depict currents 
circulating around the vortex core. The upper part shows the 
torque as a function of the line charge q(T, H) 



trolled by the Lorentz force f = (</>o/c)[J x t] exerted 
by the current density J per unit vortex length where 
t(z) = d s Y/\d s r\ is a tangent unit vector along the vor- 
tex helix parameterized by r = (u cos Qz, u sin Qz, z) and 
ds = dzy/l + Q 2 u 2 . Transport current distorts the he- 
lix, yet the net Lorentz force F = (<po/c) J Q L [J x t]ds = 
(4>qL/c)[J x z] is independent of chirality. The Lorentz 
forces acting on a helical vortex also produce the torque 
t = Jq[y x f]ds/L absent for a straight vortex. Substi- 
tuting here f = (0o/c)[J x t], we obtain that the uniform 
current density J_l perpendicular to the helix produces 
the net torque per unit length along the z— axis 



r =(f> u 2 [Jx x Q]/2c, 



(15) 



as illustrated in Fig. 3. The net torque K = tBV/4>o ex- 
erted by closed magnetization current loops vanishes, but 
a uniform current I flowing along a film strip of length 
L in a perpendicular magnetic field results in the global 
torque directed along the y-axis: 



K c = -u 2 QBIL/2c 



(16) 



Here K c (g) exhibits the behavior characteristic of the 
second order phase transition: K c = if q < q c and 
K c oc 1 — {q c /<l) 2 for q > q c even for the field H di- 
rected along the symmetry axis (see Fig. 3). This distin- 
guishes K c from the conventional torque K a = [M x H] 
of tilted straight vortices in a uniaxial superconductor 
for which K a vanishes at H||c. To estimate the mag- 
nitude of K c , we compare it with K a (ff) — VH(f>o(l — 
7 - 2 ) am20ln[r 1 H c2 /Hee}/6ATr 2 X 2 ee for H inclined by the 
angle 8 relative to the c-axis where eg = (cos 2 6+ 
7~ 2 sin 2 #) 1 / 2 [19]. For uQ ~ < > 1, we obtain that 
K c exceeds K a at J ~ Jo^tf/S, < Jo f° r any 9 where 
Jo = c</> /167T 2 A 2 £ is of the order of the depairing cur- 
rent density. Thus, the sensitive torque magnetometry 
could be used to detect twisted vortex structures. 

Helical vortices for H inclined with respect to the c axis 
may interfere with the chain and kinked vortex structures 



in layered superconductors in tilted magnetic fields [20|. 
Twisted vortex state may also affect the spiral instabil- 
ity caused by longitudinal currents in the Lorentz force 
free configurations J 1 1 H [l2j , resulting in asymmetry of 
the c-axis critical currents parallel and antiparallel to the 
twist pitch. One could also expect manifestations of the 
helical overdamped soft modes at q ~ q c in the Joseph- 
son plasma resonance in layered superconductors at H||c, 
and the effect of the chiral mixed state on electrodynam- 
ics and the magneto-optical Kerr effect 2l| . 

In conclusion, vortex charge can result in helical vortex 
instability which can enforce a spontaneous macroscopic 
twist of the vortex lattice. This can manifest itself in 
electrodynamic and thermodynamic properties of layered 
superconductors . 
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